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We study responses to metric perturbation in topological insulator models. In this paper we 
introduce a novel quantity, Hall viscosity to particle density ratio, which is analogous to the viscosity 
to entropy ratio suggested by AdS/CFT correspondence. This quantity corresponds to the filling 
fraction which characterizes the quantum Hall states, and hence it could be discussed in Haldane's 
zero-field quantum Hall system. We also consider dissipationless viscosity in the time reversal 
invariant quantum spin Hall system, which is related to a spin current. 



I. INTRODUCTION 

Recent interest in topological aspects of electron states 
has led to a large number of experimental and theoretical 
works. The study of the quantum Hall effect (QHE) [1] 
has developed the idea of the adiabatic geometric phase 
to understand their topological properties. The most re- 
markable one is that it is dissipationless transport be- 
cause a dissipative longitudinal current vanishes in the 
QH phase. However such a dissipationless transport is 
not only the Hall current. Actually, in a time reversal 
broken system, a dissipationless viscosity coefficient re- 
lated to transverse momentum transport can be discussed 

0. 0| while the conductivity corresponds to charge trans- 
port. It is called "asymmetric viscosity" [3j, "Hall vis- 
cosity" [1, Q , and also "Lorentz shear modulus" d, HI? 

1. To obtain the Hall viscosity, we can use the Kubo 
formula with a correlation function of a stress tensor as 
well as the dissipative shear viscosity. It means that the 
viscosity indicates a response to a metric perturbation. 
On the other hand, as the case of the Hall conductivity, 
the Hall viscosity is also given by the Berry curvature 
on the parameter space. For the viscosity, it turns out 
to be the moduli space of the torus 0], while we should 
consider the Brillouin zone for the conductivity [t| . It re- 
flects that the metric of the system can be parametrized 
by the modulus of the torus. 

In earlier works, mainly the Hall viscosity of the low- 
est Landau level (LLL) has been computed for non- 
interacting Q and interacting electron states, such as the 
Laughlin state and generalized conformal blocks [3, ID, Q . 
They correspond to the integer QHE (IQHE) and the 
fractional QHE (FQHE), respectively. The adiabatic 
property for the higher LL state has been actually consid- 
ered in [lfjj, but the connection with the Hall viscosity 
has not been explicitly referred in that paper. All the 
previous results are restricted to the non-relativistic sys- 
tems although there are various applications of the rela- 
tivistic Dirac fermion in the context of condensed-matter 
physics. Thus in this paper we will discuss the Hall vis- 
cosity of the Dirac fermion, and consider applications to 
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graphene [TT| and topological insulators [l2| , which have 
been remarkable topics in recent years. 

To discuss topological insulator models, we will start 
with Haldane's zero-field QH model [liij ]. In the topo- 
logical phase of the Haldane model, which turns out to 
be v = ±1 QH phase, density of states itself vanishes 
because no net external magnetic field is applied to this 
model. Indeed incompressible electron fluid cannot be re- 
alized since there are no Landau levels. As a result, the 
Hall viscosity also goes to zero even in the topological 
phase because viscosity is the hydrodynamical quantity, 
and proportional to the density or the magnetic field. 
It means that the naive definition of the Hall viscosity 
cannot capture the topological property of the Haldane 
model. Thus we have to consider a novel viscosity as- 
sociated with one-particle state, which can characterize 
the topological phase of the Haldane model, analogous 
to the filling fraction. It will be shown that a hint to re- 
solve this difficulty can be found in string theory. Recent 
progress of string theory has provided us with a remark- 
able method to approach the strongly correlated system, 
which is called AdS/CFT correspondence [l4j . Indeed 
there are a large number of applications ranging from 
particle physics to condensed-matter physics. Then the 
shear viscosity for the strongly coupled system was ac- 
tually calculated in (l5j . and universal behavior of the 
viscosity to entropy density ratio, rj/s = fi./(47rfce), was 
proposed in . It suggests that a ratio of viscosity is 
more fundamental quantity than viscosity itself, thus we 
will consider a ratio for the Hall viscosity and discuss 
some analogous properties. 

Recently it has been shown that the topological phase 
can be observed even in the time reversal invariant sys- 
tems [l7| , which is characterized by Z<i invariant [HI . Its 
discovery has a significant influence on both of experi- 
mental and theoretical sides of condensed-matter physics. 
The effective theory of such a system is well described by 
the Dirac fermion, and thus we will also generalize the 
discussion of the Hall viscosity to the time reversal in- 
variant systems. However, as discussed in Q, the Hall 
viscosity is observed only in the time reversal broken sys- 
tem. Thus we will consider another kind of viscosity in 
such a system, which is analogous to the spin Hall con- 
ductivity. We will show it characterizes a stress induced 
by a spin current while the Hall viscosity gives the rela- 
tion for a current of particles. 



2 



This paper is organized in the following way. We be- 
gin in section im with the theory of Dirac fermion as pre- 
liminaries to discuss topological insulators, and then its 
application to graphene system is discussed. In section 
IIIII we treat topological insulator models with/without 
time reversal symmetry. We introduce a novel quantity, 
a ratio for the Hall viscosity, to discuss these models, and 
some common properties with a similar quantity, which 
is considered in the context of AdS/CFT correspondence. 
We also discuss a spin Hall analogue of the viscosity ra- 
tio in the time reversal invariant systems. Section IIVI is 
devoted to summary and discussions. 



metric, is required for computing the Hall viscosity of the 
unperturbed state. This means that the Hall viscosity is 
not given by the integral of the curvature, and is not 
expected to be quantized. 



A. Dirac fermion 

As a preliminary to discuss the topological insulator 
models, we now consider the Hall viscosity of the two 
dimensional Dirac fermion model with external magnetic 
field. 

The unperturbed Dirac Hamiltonian is given by 



II. HALL VISCOSITY AND DIRAC FERMION 

Let us start with how to calculate the Hall viscosity. 
As the case of the usual dissipative viscosity, the Hall 
viscosity is obtained by response to metric perturbation, 
which is related to the shear stress mode. We now intro- 
duce the metric perturbation as follows: 



J 

9ij = — 

T2 



1 Tl 

n \r\ 2 



1 



-Tl 



Jt 2 V ~ T i 1 



(1) 

Here r = n + ir 2 is a complex form of the perturbation 
parameter, which is interpreted as the modulus of the 
torus, and g 13 stands for the inverse of gry. The Eu- 
clidean metric is recovered with t\ = 0, t 2 = 1 and 
J = 1. Then the stress tensor is given by response to 
the metric perturbation, Ty = —2(dH/dgij). Accord- 
ing to the Kubo formula, one obtains the usual viscosity 
from a correlation function of ([T xy , T xy \) . Furthermore, 
in the presence of external magnetic field, the time rever- 
sal symmetry of the system is broken, and thus another 
transport coefficient can be considered Q , which we call 
Hall viscosity. 

The Hall viscosity is actually given by correlators, such 
as ([T xx , T xy \). However it is shown in Q that the quan- 
tity is related to the Berry curvature on the moduli space 
of the metric perturbation parameters (ti,T2), while the 
Hall conductivity is given by that on the momentum 
space. The Hall viscosity is given by 



^ = 8 Im 



dg a 



dg x y 



= 2 Im 



dn 



0t 2 



.(2) 

To recover the physical units, it is required to multiply 
the factor of h/L 2 . This is just a contribution of an 
eigenstate \^> a ) to the viscosity, where a denotes a label 
of the particle state, e.g. the Landau level index. In 
general, we should consider a whole contribution of all 
the cigenstaes under the Fermi surface. 

Here we remark the difference between the Hall viscos- 
ity and conductivity. The Hall conductivity is quantized 
because it is given by the integral over the Brillouin zone 
of the Berry curvature On the other hand, the curva- 
ture on the specific point, corresponding to the Euclidean 



H = V {p x <7 x +Py(Jy) . 



(3) 



where v is the Fermi velocity. To consider the modifica- 
tion of © under the metric perturbation (JT|), we should 
find a representation of the Clifford algebra satisfying 
{<5i,CTj} = 2gij. 

The metric ([TJ can be written with the "zweibein" , 
9ij = e t el jdab, where S a b is the Euclidean metric. In the 
matrix form, it becomes 
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This helps us to obtain the representation of the Clifford 
algebra: they are given by a linear combination of the 
Pauli matrices, 



J / 1 

r 2 V r i T 2 



(5) 



Hence the Dirac Hamiltonian with the metric ((T|) yields 



n 







-irp x 



lTp x - 
iPy 



>Py 



(6) 

This is just the Dirac operator in a flat, but non-standard 
metric spacetime (see e.g. [l9l|). 

To discuss the Landau quantization for this model, let 
us take into account the magnetic field. We choose the 
Landau gauge A = (—By,0), and thus the momentum 
is modified as pj = —idj + Aj. Then the Hamiltonian 
becomes 

H = i 



a 
at 



(7) 



with uj = v\/2B / J . Here the creation and annihilation 
operators, satisfying [a, a^] = 1, are defined as 



1 



1 



V2^B 



(irp x - ip y ) , a f = 7== (~irp x + ip y ) 

(8) 

We then search for eigenstates of the Hamiltonian ^ . 
The Dirac zero mode, corresponding to the LLL state, 
is obtained by the equation of acj)^ = 0. Here a labels 
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the degenerated LLL state, a — 1,---,N, because the 
magnetic field is quantized as J F = 2irN. Then the 
higher LL states are also given by 



V n\ 



(9) 



Hence the eigenvalues and the eigenstates of the Hamil- 
tonian yield 



e (±n) =±w\/S for n = 0,1,2, •• 



$(°) = 



do) 



for n = 0, 



(±n) 




for 



1,2, 



(10) 
(11) 

(12) 



where the normalization condition for the LL state is 



given by 



dn) 



= 5 a p. As discussed in [T(j, the 



LLL wavefunction can be written in a form of 



,(0) 



— 1 /2 

M a ipa where ip a is the theta function depending on 
either of r or f. The property of the theta function yields 
the normalization constant M a = (t 2 )^ 1 ^ 2 . In general, 
the Berry curvature on the moduli space is related to the 
holomorphic structure of the eigenstate. Actually the 
LLL state consists of the (anti-) holomorphic and non- 
holomorphic parts. The latter one is just the normaliza- 
tion factor depending on both of r and f. As a result, 
it is shown that the contribution to the curvature comes 
only from the normalization factor [H, E| H| . Actually the 
zero mode, in this case, turns out to be anti-holomorphic 
function of f. According to the Berry curvature of the 
LLs obtained in [l(| , we obtain that of the eigenstates of 
0, 



r. 



(0) 
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(13) 



(14) 



Therefore the Hall viscosity is obtained from the curva- 
ture on the moduli space with the degeneracy of each LL 
state. The LL state is TV-fold degenerated, but the zero 
mode turns out to be half filled because of the parity 
anomaly (2(j. Thus contributions to the Hall viscosity 
coming from the curvature (|14[) is given by 



(™) 



hp/8 (n = 0) 
\n\hp/2 (n = ±1,±2, 



(15) 



B. Application to graphene 

The most desirable application of the relativistic (2 + 
l)-dimensional Dirac fermion theory in condensed- matter 
systems might be the graphene system [ll|, H3|. The 
low energy behavior of graphene is well described by two 
Dirac fermions at two valleys, K and K' points, and each 
valley has up and down spin degrees of freedom. Hence 
we often treat it as the four flavor Dirac fermion system. 
Indeed the anomalous QHE of graphene is well explained 
by the four flavor theory [211 ] . 

We shall discuss the Hall viscosity of the graphene sys- 
tem. First of all, it should be noted, as discussed in [j, 
it might be nonsense for the system where translational 
symmetry is violated at short length scales, since vis- 
cosity is the transport coefficient related to momentum 
transport. Therefore we have to consider the system such 
that it is well described by the continuum effective the- 
ory, i.e. the magnetic length, which characterizes the 
extent of the wavefunction, is much larger than the lat- 
tice spacing of the system. Furthermore, in the case of 
graphene, we also deal with the low energy region where 
the Dirac fermion description is suitable. 

We then introduce the effective Hamiltonian for the 
monolayer graphene [20| . 



T~Ln\ — w n 



7T 
7ft 



(16) 



where the cyclotron frequency for this system is de- 
fined as uj m = -\/2(v / 3/2)ai7As/(fi^_B), and jab is the 
hopping amplitude between the nearest neighbor sites. 
Here the annihilation/creation operators are defined as 
7r = (i£fp x - ip y )/y/2T 2 B, 7r^ = (—i£rp x + ip y )/y / 2T 2 B 
satisfying [tt, tt'] = £ where £ = 1 and £ = — 1 refer 
to the K and K' valleys, respectively. This is just the 
model discussed in section III A[ and we can apply the 
result of the Dirac fermion to this system by considering 
the metric perturbation. However, since the Dirac zero 
mode is shared by both of K and K 1 valleys reflecting the 
anomaly, we must not doubly count the contribution of 
the zero mode. Hence the Hall viscosity of the graphene 
system yields t? h = ±4(|n|(|n| + l)/4+ 1/8) ftp for n = 
0,±1,±2,---. The sign corresponds to whether it is as- 
sociated with particle states or hole states as the anoma- 
lous Hall conductivity of the graphene [2l| . On the other 
hand, in the non-interacting and non-relativistic system, 
it becomes ?7h = n 2 hp/4 for n — 1, 2, • • •, since the con- 
tribution of each LL is give n by tj^ = (n+ l/2)hp/2 for 
n = 0, 1, • • •, as shown in 



ven Dy i 
[4, 10]. 



where p = N/L 2 = l/(2n£ 2 B ) is the density of particles HI. HALL AND SPIN HALL VISCOSITY RATIO 

and the magnetic length is defined by £b = \/h/(eB). 

We remark the result (TT51) depends on the level index n, In this section we then apply the argument of the Hall 

while the contribution to the Hall conductivity does not. viscosity to some topological insulators [12]. All the topo- 

It has been shown this level dependence indeed repro- logical insulator models have a band gap, but there are 

duces the classical result = k-^Tpl (2w c ), where uj c is some ways to generate a gap: it is generated by mag- 

the cyclotron frequency @,13]- netic field in a time reversal symmetry broken system, 
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and spin-orbit interaction causes it in a time reversal 
symmetric system. Although the origin of the gap is dif- 
ferent, indeed topological aspects of such systems are well 
described by the effective Dirac fermion model discussed 
in section [TT1 

A. Hall viscosity ratio in topological phase 

First we consider the model proposed by Haldane 13:]. 
The remarkable property of this model is that the topo- 
logical phase is caused by only breaking time reversal 
symmetry while there is no net external magnetic flux. 
Indeed it is shown that anomalous contribution of the 
Dirac zero mode, relating to the parity anomaly, plays 
an important role in it. 

Now let us discuss how to obtain the Hall conductivity 
for the model without external magnetic field. It is given 
by a derivative of charge density with respect to mag- 
netic flux, <7h = dp/dB, which is called Widom-Stfeda 
formula [22], HH and available even in the zero field limit. 
It can be interpreted as l'Hopital's rule to evaluate an 
indeterminate form of the filling fraction in the limit of 
B4 0. Indeed the charge density becomes p = ±e 2 B/h 
by taking into account infinitesimal magnetic field in a 
certain phase of the Haldane model. Therefore it turns 
out to be v — ±1 QH phase. 

We then consider the Hall viscosity of the topological 
phase in the Haldane model. According to the result of 
the Dirac fermion model discussed in section [Hi it yields 
T^H = i (1/4) Zip in the topological phase with infinitesi- 
mal magnetic flux. We now remark that not charge, but 
particle number density is used in the expression of the 
Hall viscosity (1151) . However, taking the zero field limit, 
the Hall viscosity goes to zero as the particle density van- 
ishes. This reflects that viscosity is the hydrodynamical 
quantity which depends on the density. In the Haldane 
model the incompressible fluid state on the LLs cannot be 
found, and thus the naive definition of the Hall viscosity 
cannot capture the topological property of the model. To 
resolve this difficulty, we now introduce a novel quantity, 
a ratio of the Hall viscosity to particle density, 

m = m/p- (17) 

This is similar quantity to kinematic viscosity, which is 
defined as the ratio of viscosity to mass density. The 
reason why we prefer the quantity (1171) to kinematic vis- 
cosity is that it is not well defined for the massless Dirac 
fermion. Hence, by using l'Hopital's rule as the case of 
the Hall conductivity in the zero field limit, which can be 
interpreted as alternative formula for the Hall viscosity 
ratio of Widom-Stfeda formula, we see that it can takes 
non-zero value, 

m = ±\- (is) 

In this sense, it stands for t)h = sh/2 with effective one- 
particle spin s because the Hall viscosity is represented 



as ?7h = 

We now comment universal behavior of the Hall vis- 
cosity for the LLL state. The Hall viscosity of the LLL 
state, which is discussed in [3, 6], and that of the Dirac 
zero mode shown in (j!5[) are given by 77H = hp/ A, yield- 
ing ?)h = h/A. Furthermore it is shown in [4j, |8j that 
the Hall viscosity of v = 1/m Laughlin state becomes 
t/h = mhp/A with p — l/(2mirl 2 B ). Although the value 
of the Hall viscosity itself is different from the above, 
after all dependence on the magnetic field provides the 
universal value, r/n/B = e/(8n). On the other hand, 
this universality is not observed in the excited state. In- 
deed the Hall viscosity for more generalized conformal 
blocks is obtained in [J], rjn = (v -1 /2 + h^)hp/2, with a 
conformal weight h^. Similarly, the n-th LL state gives 

7 ?h"' = ( n + 1/2) ftp/2 as discussed before. These results 
give a larger value of the ratio r/a/B > e/ (87r) compared 
with the LLL state because an excited state has higher 
spin. 

We also remark similarity of the Hall viscosity ratio r)jj 
to the viscosity to entropy density ratio, rj/s. According 
to AdS/CFT correspondence, the usual dissipative vis- 
cosity for the strongly coupled system is obtained in [l5j , 
and its universal behavior, rj/s = h/{Ai:k^,), called KSS 
bound, is discussed in QJ]. In the context of AdS/CFT, 
the Hawking temperature of the black hole is interpreted 
as that of the dual field theory at boundary. Although 
the temperature goes to zero in the extremal black hole 
limit, the residual entropy is observed in some cases. Ac- 
tually it is discussed in [24| [H[ that this universality is 
still available in the extremal case, corresponding to the 
zero temperature limit although the entropy also goes 
to zero in the limit. On the other hand, it is shown in 
26] that the KSS bound can be violated in generalized 
holographic model, e.g. Gauss-Bonnet gravity, which in- 
cludes higher order curvature terms. In the case of Hall 
viscosity, the shift of the Hall viscosity, coming from the 
effect of the geometry, is also pointed out [J|. These 
results suggest the geometrical effect should not be ne- 
glected for arguments on viscosity. 

As discussed above, there are some common properties 
in both of the viscosity ratios. However these universal 
behaviors mean that viscosity itself might not be uni- 
versal quantities. The reason why is the natural unit of 
the viscosity should be H times density as discussed in 
[H . This density dependence suggests viscosity relates to 
the characteristic scale of the system, e.g. the magnetic 
length. Thus it is preferred that we consider the ratio of 
quantities rather than viscosity itself. 

B. Spin Hall viscosity with time reversal symmetry 

It is well known that the Hall conductivity occurs if 
and only if time reversal symmetry is broken. As the 
case of conductivity, it is pointed out in Q that the Hall 
viscosity can be also observed in a system without time 
reversal symmetry. On the other hand, spin Hall conduc- 
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FIG. 1: Stress induced by (a) chiral current in QH state and 
(b) helical current in QSH state, which is equivalent to the 
spin current. Solid (dotted) lines represent currents of up 
(down) spin. 

tivity induced by spin-orbit interaction is found in time 
reversal symmetric systems [13, H^. We then consider 
application of the idea of the Hall viscosity to the time 
reversal system, and spin Hall analogue of viscosity. 

The most fundamental model for the topological in- 
sulator with time reversal symmetry is the Kane-Mele 
model, which is proposed in [Tt], Gil . In this model the 
spin-orbit interaction is taken into account on honeycomb 
lattice, which corresponds to the next nearest neighbor 
(NNN) hopping term of the Haldane model. Thus the 
Kane-Mele model can be interpreted as the doubled Hal- 
dane model, because the signs of the NNN hopping am- 
plitude is opposite for up and down spins. Actually the 
time reversal symmetry of the Haldane model is explic- 
itly broken, which is equivalent to up/down spin sector 
of the Kane-Mele model. However the Kane-Mele model 
preserves the time reversal symmetry by considering both 
of up and down spin sectors. Hence a quantized Hall con- 
ductivity of the Haldane model induces a spin current, 
J s = (h/ (2e))(Jj- — Jj.), thus it yields a spin Hall conduc- 
tivity <7 s h = e/(27r). In this sense, we can consider spin 
Hall analogue of the Hall viscosity ratio 

^H = i(4 t) -# ) )=^ ^ 

which we will call spin Hall viscosity ratio. As the case of 
the Haldane model, the quantity (|19l) is just for the one- 
particle state. Although some generalizations of topolog- 
ical insulators for correlated electron systems have been 
considered [29|, [3(| , an effect of many-body interaction is 
not yet obvious. If certain aspects of highly correlated 
topological insulators are understood, we would be able 
to discuss the spin Hall viscosity for such systems. 

We then study a physical meaning of the spin Hall vis- 
cosity. In the presence of magnetic field and a current 
the Hall viscosity gives a relation between the current 
and the stress tensor, T xx = —T yy = rmd y v x , where v x 
stands for the current in x-direction Q. In the case of 
the spin Hall system, it corresponds to T xx — —T yy = 

Vsn(d y v x ^ - d y Vx) = r/ s nd y v s/x where v s , x represents 
x-component of a spin current. This means the spin cur- 
rent induces the stress in the Kane-Mele model as shown 
in Fig. [T] while the charge current does not. It is dissipa- 



tionless because the stress is perpendicular to the current 
as the cases of Hall and spin Hall conductivity. 

IV. DISCUSSIONS 

In this paper we have discussed the Hall and spin Hall 
viscosity in some topological insulator models. In section 
III1 the metric perturbation for the Dirac fermion model 
is studied. The Berry curvature on the moduli space, 
which gives the Hall viscosity of the Dirac fermion, is 
computed, and its application to the graphene system 
is also discussed. In section |HT] we have proposed the 
novel quantity, the Hall viscosity to particle density ratio, 
to study the viscosity of Haldane's zero-field model for 
QHE. It is analogous to the ratio of viscosity to entropy 
density, which is discussed in the context of AdS/CFT 
correspondence, and some common properties are found. 
We have also studied the spin Hall viscosity of the Kane- 
Mele model. Although the Hall viscosity is observed only 
in the time reversal broken systems, the spin Hall viscos- 
ity can be found even in the time reversal invariant sys- 
tems. It is then pointed out that the spin current induces 
the dissipationless stress in the topological phase of the 
Kane-Mele model. 

We now comment some open problems on this topic. 
One of them is the relation between the Hall viscosity and 
the edge state. It is well known that the edge state plays 
an important roll on the transport phenomena in the 
QH and QSH systems: the Hall conductivity is actually 
obtained by different ways, due to the bulk/edge corre- 
spondence [31| . In this sense, the Hall viscosity should be 
computed from the edge state. Next is a generalization 
to the three dimensional topological insulators, which are 
also characterized by Z 2 invariant [32ll34l ]. e.g. the Fu- 
Kane-Mele model [32| . Three dimensional anomalous vis- 
cosity induced by the external field was indeed considered 
in Q , and its classical computation was also given. Since 
the surface state of the three dimensional model can be 
described by Dirac fermions, the results discussed in this 
paper would be related to it. Furthermore the four di- 
mensional generalization would be also discussed. The 
anomalous contribution in the four dimensional hydro- 
dynamics is recently considered in (35| . It is related to 
the triangle anomaly, and corresponds to the non-linear 
response coefficient. Indeed, in the four dimensional QH 
system (36[, the non- linear coefficient is quantized while 
the linear part is not. Therefore it is expected that the 
four dimensional anomalous contribution can be found in 
the four dimensional QH system. 
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